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ABSTRACT 



The properties of two-state nearest-neighbour cellular 
automata (CA) that are capable of density classification 
are discussed. It is shown that these CA actually con- 
serve the total density, rather than merely classifying it. 
This is also the criterion for any CA simulation of DFT. 
The effect of boundary and periodicity conditions upon 
the evolution of such CA are elaborated by considering 
linear and cyclic lattices and boundary conditions consis- 
tent with the conservation criterion. In a bounded linear 
lattice, it is possible to achieve a configuration with a sin- 
gle 01 and no 10 domain wall, or vice versa, but this is 
not possible in a cyclic lattice, where these two domain 
walls have to appear in pairs. This determines the final 
stable state of the automaton. 
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It has recently been shown by Capcarrere, Sipper and 
Tomassini 2 that certain two- state r=l Cellular Automata 
(CA), namely rules 184 and 226 in Wolfram's classification 3 , 
can perfectly classify the density, in contrast to a result 
proved earlier by Land and Belew 4 . This is due to a re- 
definition of the classification criterion for the final con- 
figuration in ref.2. The final configurations for these CA 
consist of a background of alternating bits, with a block of 
one or more adjacent O's or l's, depending upon whether 
the density of l's in the initial configuration was less than 
or greater than 0.5. No other CA performs this density 
classification with comparable accuracy: rules 57 and 99 
can classify the density with only 60% accuracy. 

In this paper I point out the properties of the CA 
that make such density classification possible and the ef- 
fect of boundary conditions upon the evolution of such 
CA. First consider the same periodic boundary condi- 
tions (sn+i = Si where Si is the value of site i) employed 
in ref.2 and define a "vacuum state" consisting purely 
of alternating O's and l's. The density classification is 
performed by the mutual annihilation of quasiparticles 
(two or more adjacent l's) and holes (two or more ad- 
jacent O's), which move in opposite directions through 
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the lattice, until there are no further quasiparticles or 
no further holes, depending upon which were fewer in 
the initial configuration. A pair of adjacent l's and a 
pair of adjacent O's (in a background of alternating bits) 
are topological solitons : one is a kink and the other an 
antikink. The topology is identical to that of quasi- 1- 
D polymers like polyacetylene, where solitons have been 
identified and simulated by CA 5 . 

During the evolution of the CA rules 184 or 226, the 
number of O's and l's are each conserved (so these CA are 
not merely classifying the density, but rather conserving 
the density). Such conservative CA were first discussed 
by Vichniac 6 . This conservative evolution is possible be- 
cause CA rules 184 and 226 can be reformulated as two- 
site density interchanges or two-site spin flips (if the bits 
are visulaized as spins). Rule 184 may be restated : 

Flip any 10 > 01; All else unchanged (1) 

while rule 226 may be restated as : 

Flip any 01 > 10; All else unchanged (2) 

This makes l's and O's (or opposite spins) propagate in 
opposite directions, unlike other conservative CA (rules 
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170 or 240) where both species drift in the same direc- 
tion. Note that rules 57 and 99 are not perfect density 
classifiers, since they donot conserve the density (they 
change a block of contiguous 0's into a block of l's and 
vice versa, at each step) . 

The final state <E> a ^ of a CA (rule 184 or 226) consist- 
ing of N sites, with average density p = 0.5, is achieved 
within N/2 steps and oscillates between a configuration 
where S{ = for all even i, = 1 for all odd i and a 
configuration where = 1 for all even i, Si = for all 
odd i (an " antiferromagnetic" state in the spin analogy). 
When the density is different from 0.5 there are solitonic 
structures, consisting of either only l's (for p > 0.5) or 
only 0's (for p < 0.5) superposed upon <E> a ^ , which prop- 
agate continuously around the lattice. 

The situation changes markedly when the boundary 
condition is no longer periodic. Another set of possible 
boundary conditions consistent with the criterion of den- 
sity conservation is : 

s = 0, s N+1 = 1 for Rule 184 (3) 

s = 1, s N+l = for Rule 226 (4) 
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With a suitably bounded linear lattice, the "antiferro- 
magnetic" state & a i t (obtained with the periodic lattice) 
is no longer stable. <E> a / t evolves further, to give a state 
consisting of a single block (<E> ) of O's and a single block 
(<E>i) of Is ("ferromagnetic" domains), which is then a sta- 
ble end point of the CA evolution. The entire evolution is 
completed within N-l steps. Writing Wj\<j){m,n) > for 
the application of Wolfram's CA rule r upon an arbitrary 
initial configuration |0(m, n) > consisting of m O's and 
n l's, for T steps, we now have 

Wl M \(j)(m, n) >= |$ (ra)$i(n) > for all T > N - 1 

(5) 

where the block of m O's appears to the left of the block 
of n l's, and 

Wg 6 |0(ra, n) >= |$i(rc)$ (ra) > for all T > N - 1 

(6) 

where the block of m O's appears to the right of the block 
of n l's. Thus VF^ 1 and W^q 1 act as ordering opera- 
tors (as in the normal ordering of creation and annihila- 
tion operators in many body theory). 

The effect of the boundary conditions upon the final 
configurations of the CA arises from the changed topology 
of the lattice. From equation (1), it is seen that the 10 
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domain wall is unstable under Wi&4 while the 01 domain 
wall is stable. In a bounded linear lattice, it is possible 
to achieve an arrangement with a single 01 and no 10 
domain wall (as in the RHS of equation 5), but this is 
not possible in a periodic lattice (which has the topology 
of a circle) where these two domain walls have to appear 
in pairs. Likewise the 10 domain wall is stable under 
W226 while the 01 domain wall is unstable. As a result 
<& a it is unstable under rules 184 and 226 in a linear lattice 
with the boundary conditions (3) or (4). 
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